Magnetotransport in a two-dimensional electron gas 
in the presence of spin-orbit interaction 

X. F. Wang and P. Vasilopoulos 

Department of Physics Concordia University 
Montreal, QC H3G 1M8, Canada 
(Dated: February 2, 2008) 

Abstract 

We evaluate the transport coefficients of a two-dimensional electron gas (2DEG) in the presence 
of a perpendicular magnetic field and of the spin-orbit interaction (SOI) described only by the 
Rashba term. The SOI mixes the spin-up and spin-down states of neighboring Landau levels into 
two new, unequally spaced energy branches. The broadened density of states, as a function of the 
energy, and the longitudinal resistivity, as a function of the magnetic field, show beating patterns 
in agreement with observations. The positions of any two successive nodes in the beating pattern 
approximately determine the strength of the Rashba term. A strong SOI results in a splitting of 
the magnetoresistance peaks and a doubling of the number of the Hall plateaus. The peak value 
in derivative of the Hall resistivity reflects the strength of the SOI. 

PACS numbers: 73.20.At; 73.20.Dx; 73.61.-r 



1 



I. INTRODUCTION 



There has been an increasing interest in zero-magnetic-field spin splitting in one- (ID) 
and two-dimensional (2D) electron systems due to the spin-orbit interaction (SOI). Such 

riQ 

systems have potential applications in spin-based transistors |2j expected to service in the 
future quantum computation. The SOI has been found also important in an unexpected 
metal-to-insulator transition in 2D |3j hole gas, in spin-resolved ballistic transport jj|, in 
Aharonov-Bohm A-B experiments 5], and in a spin-galvanic effect [6]. The analysis of the 
Shubnikov-de Haas (SdH) oscillations in magnetoresistance measurements has become the 
main method of measuring the SOI strength in such systems. 

Decades ago theoretical studies jl, |l| in 3D semiconductors found that the spin degen- 
eracy should be lifted in inversely asymmetric crystals due to the internal crystal field. 
Later, magnetotransport and cyclotron resonance measurements in a 2D hole system, in 
a modulation-doped GaAs/AlGaAs heteroj unction, showed [l| evidence of zero-magnetic- 
field spin splitting for carriers with finite momentum. Similar experiments on 2D electron 
gases, formed in a GaAs/AlGaAs inversion layer, led to similar conclusions Q]. The first 
explanation was proposed by Bychkov and Rashba [ljj employing the Rashba spin-orbit 
Hamiltonian, where the spin of finite-momentum electrons feels a magnetic field perpendic- 
ular to the electron momentum in the inversion plane. Though nonparabolicity of the bulk 
band structure of GaAs/AlGaAs could also explain the previous experimental results and 
bulk inversion- asymmetry induced spin splitting, at B=0, could dominate in heterostruc- 
tures of wide-gap semiconductors, the Rashba SOI has been considered the most appropriate 
reason for the observation of the zero-field spi n splitting in low-dimensional electron systems, 



especially in narrow-gap semiconductors 



12|. Later. Luo et al. Q investigated the SdH 



oscillations in a series of GaSb/InAs quantum wells and concluded that the lifting of the spin 
degeneracy results from the inversion asymmetry of the structure which invokes an electric 
field perpendicular to the layer. Using the Rashba SOI, they fit the experimental results 
and determined the Rashba parameter a, which describes the strength of the SOI. At the 
same time they concluded that contributions to the SOI from the bulk ~ k 3 term due to a 
crystal inversion asymmetry are of minor importance. 

Generally, the contributions to the spin splitting in the conduction band of asymmet- 
ric heterostructures result from the bulk ~ k 3 term due to a crystal inversion asymmetry 
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and from the Rashba ~ k asymmetry. Due to their different momentum dependence, the 
former dominates in wide-gap structures with small thickness whereas the later dominates 
in narrow-gap structures. It was shown ^| that the k 3 term leads to anomalous beating 
patterns while the Rashba term leads to the regular beating patterns in magneto oscilla- 
tions. Recently, the well-developed shaping technique in nanostructures has been used to 
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20, 



2l|, and 



control the SOI strength in 2D systems of different materials 
principally the SdH oscillations are used to measure the Rashba parameter n [22]. However, 
to our knowledge there are no detailed theoretical treatments of the influence of the SOI 
on magnetotransport in 2D systems. We therefore aim at developing a more realistic model 
to describe theoretically magnetotransport in systems with SOI, in which the Rashba term 
dominates, and determine more accurately the parameter a. 

In Sec. II we present the energy spectrum and the density of states (DOS). In Sec. Ill we 
present the results for the transport coefficients and in Sec. IV concluding remarks. Some 
auxiliary results are found in the appendix. 



II. EIGENVECTORS, EIGENVALUES, AND DENSITY OF STATES 

We consider a 2DES in the (x — y) plane and a magnetic field along the z direction. In 
the Landau gauge A = (—By, 0, 0) the one-electron Hamiltonian including the Rashba term 
reads 

H = ^t^ 2 + % W x (P + eA)l, + 9PbB<t„ (1) 
2m* n 

where p is the momentum operator of the electrons, m* is the effective electron mass, g the 
Zeeman factor, fiB the Bohr magneton, o = (o~ x ,o~ y ,o~ z ) the Pauli spin matrix, and a the 
strength of the SOI. 

Using the Landau wavefunctions without SOI as a basis, we can express the new eigen- 
function in the form (k x commutes with the Hamiltonian (1)) 



**» = ^^^(2/-2/ c )C»/v^ 



e*- x 52Mv-Vc) [ °\ ] /V 7 ^, n = 0,l,2,---. (2) 

C„ 
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Here <p n (y — y c ) = e~^ v ~ v ^ ^ 2l ^H n {{y — y c )/l c )/ a/ y / 7r2 n n!/ c is the usual harmonic oscillator 
function, uj c = eB/m* the cyclotron frequency, l c = (h/m*uj ( ) l l 2 the radius of the cyclotron 
orbit centered at y c = l 2 k x , n the Landau- level index, and \a) the electron spin written as 
the column vector |cr) = ^ * J if it's pointing up and ^ ° j if it's pointing down. Substituting 
Eq. (J2J) in the Schrodinger equation H^/ = E^/, multiplying both sides by (pi(y — y c ), and 
integrating over y we obtain the following system of equations (E± = ±gfiBB — E) 



i{a/l c )V2lCf_ x + [(/ + l/2)hw c + E_)Cf = 



I 



0,1,2,.-. 



(3) 



^ [(/ + l/2)tkj e + E + ]C+ - i(a/l c )y/2(l + l)Cr +1 = ^ 
This infinite system of equations is solved exactly after decomposing it into independent 
one- or two-dimensional secular equations. Denoting the new subband index by s we obtain 



[l/2hu) + E-]C~ = 0, s = 0; 



(4) 



[s-l/2)huj c + E + -i(a/l c )y/2s 
i{a/l c )V2s (s + l/2)fkv c + E_ 



Cr 



0, S =l,2,3,-... 



(5) 



Corresponding to s = 0, there is one level, the same as the lowest Landau level without SOI, 
with energy Eq — Eq — l/2ftw c — gfisB and wave function \I/q (A^) = (e lkxX ' / ' yfL^)4>o{y — 
y c ) ^ ° j . Corresponding to s = 1, 2, 3, ■ • • , we find two branches of levels with energies 



Ef = sfku c ± [El + 2sa 2 /l 2 c } 1 ' 2 . 
The + branch is described by the wave function 



(6) 



y/L x A s 



-W s (j) s ^i(y - y c ) 

<Ps(y - y c ) 



(7) 



and the — one by 



y/L x A s 



<t> s -i{y-y c ) 

-ffi s (f> s (y - y c ) 



(8) 
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where A s = 1 + V 2 S and 

V = y/Taallc (9) 

S E + ^Ei + 2sayP c ' 

The density of states (DOS) is defined by D(E) = J2 s k x a ^(-^ — )■ Assuming a Gaussian 

broadening of width T we obtain 

o p -(£-_B^) 2 /2r 2 

In Fig. 1 (a) we plot the level energies E+ and i£~ as functions of the level index s. 
For the case studied here we have E± ~ Eq . Because the level spacing of the + branch 
is larger than that of the — branch, the level energy of the + branch increases faster and 
the line through the triangles (not shown) has a slope larger than that through the circles 
(not shown) in Fig. 1 (a). Here we also notice that ~ (E£ + E£)/2, E± 5 ~ E± 3 , 
E^ 5 — {E22 + -S23V2. This difference in level spacing results directly in the modulation of 
the density of states as shown in Fig. 1 (b) and (c). As Fig. 1 (b) shows, where the level 
broadening is small, Y = 0.1 meV, the DOS as a function of the energy shows peaks of the 
same height except when levels of different branches have the same value and higher DOS 
peaks appear. For wider level broadening, as shown in Fig. 1 (c), with Y = 0.5 meV, the 
DOS is modulated and shows a beating pattern. The nodes of this pattern appear when a 
— branch level is located near the middle between two + branch levels; thus, the first node 
appears near the E^ level and the second node near the E^ one. The maximum oscillation 
amplitude appears when two levels of different branches are degenerate, e.g., at E± 5 ~ E± 3 
here. 



III. TRANSPORT COEFFICIENTS 
A. Analytic results 

For weak electric fields E U) i.e., for linear responses, and weak scattering potentials the 
expressions for the direct current (dc) conductivity tensor a^ u , in the one-electron approxi- 
mation, reviewed in Ref. j^, reads a^ u = a A + with \i,v = x, y, z. The terms <j a and 
cx™^ stem from the diagonal and nondiagonal part of the density operator p, respectively, in 
a given basis and (J M ) = Tr(pJ M ) = cr flu E u . In general, we have a A = a^J + a c °^. The term 
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B(T) 

FIG. 1: (a) Subband energy E s versus index s. The triangles are for the + branch and the circles 
for the — branch, (b) Energy (right scale) versus DOS with a subband broadening T = 0.1 meV. 
(c) The same as in (b) but with T = 0.5 meV. The other parameters are g = 2, B = 1 tesla, 
m* = 0.05, and a = 10- n eVm. 

a^J describes the diffusive motion of electrons and the term a c °^ collision contributions or 
hopping. The former is given by 

< f = p 4- E - mmE^, (ii) 

where ( = (s,a,k x ) denotes the quantum numbers, i>£ = (Cl^vIC) i s the diagonal element 
of the velocity operator and f(e) the Fermi-Dirac function. Further, t^(E^) is the 
relaxation time for elastic scattering, (3 = 1/ksT, and So is the area of the system. The 
term of,?' can be written in the form 



a 2 poo poo 

< = VE / d£ / ^[ £ -^(^)]5[ £ '-^;(^)]/(£)[i-/(£0]Wcc(^^)(yc-ycO 2 , 

C C' J —oo 

(12) 

where = (CMC); W^r(e,E') is the transition rate. For elastic scattering by dilute impuri- 
ties, of density Nj, we have 

W cc (e,e>) = ^El^l^CC'WI 2 ^-^)^,^, (13) 
q 
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where u = l 2 q 2 /2 and q 2 = q 2 + q 2 . U(q) is the Fourier transform of the screened impurity 
potential U(r) = (e 2 /47re e)e _fcsr /r, where e is the static dielectric constant, eo the dielectric 
permittivity, and k s the screening wave vector. 



(14) 



2e e {2u/l 2 + k 2 y/ 2 

In the situation studied here the diffusion contribution given by Eq. (fTTj) vanishes because 
the diagonal elements of the velocity operator vanish. Neglecting Landau-level mixing, i. 
e., taking s' = s, and noting that o c °/. = a™ 1 , = (So/ 2%) J °° qdq = (So/2nl 2 ) J °° du, and 
Efe = (So/2tt1 2 ), we obtain 

.coz Njf3e 2 



a 



mi 



2M 



i>oo poo 2 

J2 du ds[S(s-E^} 2 f(e)[l-f(s)}u(^M) |F S »| 2 M , (15) 



where 



\F7.(u) 



F+(u) 



{L s ^(u)+V 2 s L s (u)Ye- u /A 
{V 2 s L s ^(u)+L s (u)} 2 e- u /A 



The exponential e " favors small values of u. Assuming b = k 2 l 2 /2 ^> u 24|, we may 
neglect the term 2u/l 2 in the denominator of Eq. (fUjl and obtain 



a 



III) 



Ntfe 2 
Aithb 



2ee 



/oo 
de[5(e-E°)ff(e)[l-f(e)\r s , 
-oo 



where 



[2s ± 1)X>* - 2s£> 2 + 2s± 1}/A 2 S 



(16) 



(17) 



The impurity density Ni determines the Landau Level broadening T = W^> (e, e') / h. 
Evaluating W^i(e,e')/h in the u — » limit without taking into account the SOI, we obtain 
Ni « 47r[(2ee /e 2 )] 2 r/^. 

The Hall conductivity a™^ is given by 



(7 



2z/Ie 2 ^ 1 - e^- E 'c ] 

" - " E / " /Wl < C I ** I C >< C I I C > /F , C * C (18) 



S 
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FIG. 2: Conductivity function of the magnetic field B. The dotted curve is for a = 

eVm and the solid one for a = 1.2 x 1CP 11 eVm. The inset shows the oscillations between the 
fourth and the fifth node. 



The evaluation of Eq. (fT%|) proceeds along the lines of Ref. |2j| using the the matrix elements 
< C I v n I C >> A 4 = x ; Z/> given in the appendix. Taking — e u n ~ ftu; c , leads readily to 



where 



a 



nd 
xy 



e" 



j> + 1) [B s (ft - + C s {f- +1 - /- 2 )] 



s=0 



(19) 



/I 2 ,4 2 



2m*aV 



s+l 



h^/s + 1 m*/ c 



e /2 , 2m*aV 8+1 

h 2 uj c \/ s + 1 



1 



el 



(20) 



(21) 



A 2 A 2 

here 6 S = 1 + ^/(s + l)] 1 / 2 ^!?^! and 0' s = 1 + [(s + 2)/(s + l)] l / 2 V s+l V s+2 . We notice that 
if a = 0, we have Q s = Q' s = A s+ i = A s+2 = 1 and Eq. (fT^j) becomes cr"^ = (e 2 /h) J2t =0 f s , 
i.e., the conductivity expression pertinent to the integer quantum Hall. [25| 

The resistivity tensor, p^ is given in terms of the conductivity tensor. We use the 
standard expressions p xx = a yy /S, p yy = cr xx /S, p yx = p xy = -a yx /S, where S = a xx a yy - 



@xy@yx • 



s 



B. Numerical results 



In the numerical evaluation of the conductivity we assume that the S functions appearing 
in Eq. are broadened and replace them with the Gaussian function (1 / '\^2nT)e~ x ' 7( 2r2 ). 
Further, if not otherwise specified, we use the following parameters: T = 0.4K, T = 1.5meV, 
g = 0, m* = 0.05, iV e = 4x 10 16 m~ 2 . In Fig. 2 we plot in the small u limit, cf. Eq. 
f!16|) . as a function of the magnetic field. The dotted curve shows a^J in the absence of SOI 
and the solid one in its presence with a = 1.2 x 10 -11 eVm. For low magnetic fields and 
weak a the conductivity decreases quickly with a and saturates around a = 2 x 10 _11 eVm. 
The typical beating pattern appears when the subband broadening is of the same order as 
the Landau-level separation. At high magnetic fields, the effect of SOI is weakened and the 
beating pattern is replaced by split conductivity peaks. The latter approaches that without 
SOI when the magnetic field becomes very strong. One of the segments, with a typical 
beating pattern between the fourth and the fifth node, is shown in the inset. From Fig. 1 
we see that the mth node is located near the nth Landau level when E~ +m ~ (E^_ 1 + E^)/2. 
Using this expression for large n and small m (n ^> 1), we can obtain the ratio of the Landau 
index over the magnetic field. The result is 

n m (2m - l)(2m + 3)eh 2 + 8gm*he - Am* 2 g 2 h 2 e 

(22) 



This leads to 



B m 32m* 2 a 2 



n m+ i n m (m+l)h 2 e 



(23) 



B m+ i B m Am* a 2 
If we keep the electron density N e constant and use the definition of the filling factor v = 
N e 2irl 2 , we can approximate B m by B m = 7rhN e /en m , and obtain 

1 1 (m + l)he 2 

(24) 



B 2 m+1 Bl 47im* 2 a 2 N t 



and 



2 (m + l)nh 3 
m+1 m ~ 4m* 2 a 2 



(25) 



Eq. (J2"2j) and (J23j) can be used to estimate the Rashba parameter a. For instance, using the 
inset of Fig. 2 provides = 71 and 715 = 87. Experimentally, the SdH oscillations in the 
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FIG. 3: The resistivity 

Pxx ^ 

function of the magnetic field B for different values of the parameter 

a as indicated. 

resistivity of a 2D system, in the presence of SOI, are usually viewed as resulting from a 



2D system with two subbands 



pre 

m 



15. 21] with the SOI splitting at the Fermi level A^ = 2akp 
serving as the subband separation. Following this line of reasoning, we can also analyze 
the results shown in the inset of Fig. 2. The SdH frequency difference between the plus and 
minus oscillations is m x B m ~ 4.8 Tesla and corresponds to a carrier density difference 
AN = 1.16 x 10^ m~ 2 . This leads [18] to a = Hk F AN/(2m*N e ) = 1.1 x lCr n eVm. 

In Fig. 3 we plot the resistivity p xx for different strengths a as a function of the magnetic 
field. With the increase of a each resistivity peak becomes lower and gradually splits into 
two peaks. However, the shape of the gaps is not affected by SOI. We also notice that all 
peaks retain almost the same form after splitting. 

Figure 4 shows the Hall resistivity p xy versus the magnetic field B. For strong magnetic 
fields we see the integer quantum Hall effect plateaus at h/ne 2 , where n is a integer. In 
the presence of SOI, one more plateau with value 2h/{2n + l)e 2 appears between every two 
plateaus of order n and n + 1. The size of this new plateau increases with a. It is worth 
noting that these extra plateaus require rather strong a and may easily shrink or disappear 
if disorder is included in the calculation of the Hall resistivity. In the lower inset of Fig. 4 
we plot the derivative dp xy /dB as a function of B. Each peak, corresponding to a sharp 
jump of the resistivity, splits into two peaks which separate from each other, by a distance 
AB, with increasing a. The dependence of AB on a is plotted in the upper inset. The split 
increases slowly for small a and saturates at about a = 2 x 10~ n eVm. 



10 




10 15 20 25 

BCD 

FIG. 4: Hall resistivity p xy as a function of the magnetic field B. The different curves correspond 
to different a and are marked as in Fig. 3. The lower inset shows the derivative of p xy with respect 
to B versus B. In the upper inset the difference AB between the values of the two peaks in this 
derivative, into which the a = peak near B = 24 tesla splits, is shown as a function of a. 

C. Comparison with the experiment 



In the following, we will analyze, using Eq. (|24|) . two typical measurements of the SdH 
oscillation in InGaAs/InAlAs heterostructures, assuming the Rashba term dominates the 
contribution of the observed zero-field spin splitting. Ref. 26J provides results for two 
samples, A and C. For sample A, with effective mass m* = 0.046 and sheet density n s = 
1.75 x 10 12 cm -2 , the positions of the first six nodes are, respectively, at fields, 5i=0.873T, 
5 2 =0.46T, 53=0.2917, 5 4 =0.227T, 55=0.1837, 56=0.1537. From the positions of any 
two successive nodes, B m and B m+1 , we extract the Rashba parameter a. The results are 
shown as full triangles in Fig. (j3J); as shown a fluctuates and converges to the average 
value a = 3.7 x 10~ 12 eVm with increasing node number m. The consistency of the a 
values extracted from different nodes convinces us that the Rashba term is the main cause 
of the beating pattern here. It may be that bulk SOI contributes also and results in the 
variation of a when different nodes are used. The calculated spin splitting at the Fermi 
level is 5p = 2akp = 2.45 meV and is the same as the extrapolated result from Ref. 
The same analysis has been done for sample C, with n s = 1.46 x 10 12 cm -2 and 5!=0.65T, 
5 2 =0.312T, and B 3 = 0.204T; the results are shown in Fig. El as squares. The zero-field 
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FIG. 5: Strength a of the Rashba term, as a function of the observed [2(| node number m, 
extracted from Eq. (|24j) and measured node positions B m and B m+ \ for sample A (triangles) and 
sample C (squares). The inset shows our calculated p vs B beating pattern for sample A. The 
dotted and solid lines are guides to the eye. 
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FIG. 6: Strength a of the Rashba term, as a function of the applied gate voltage, extracted 
from Eq. (|24|) (squares) and pertinent to the results of Ref. The latter are shown by the 

open (solid) circles when the first (second) nodes are fitted as in Ref. The inset shows 

our calculated p vs B beating pattern for V g = 0.3V, T = 0.4K. The dotted curve, produced by 
a = 7.04 - 2.26Vg + 0.87F 2 , is a fit to our results (squares). 



spin splitting at the Fermi level is 1.7meV and is close to the value 1.5meV given in Ref. 



26]. In the inset of Fig. (jSJ) we show the calculated diagonal resistivity for sample A of Ref. 

261 ] at temperature T = 0.5 K. A magnetic-field-dependent subband broadening is adopted, 
T = Tq\/B, with r = 0.68 meV/T 1 / 2 . The second node is well fitted, while the first node 
appears at a slightly higher magnetic field and 40 oscillations are enclosed between them 
whereas the number observed in Ref. |26j is 35. 



12 



In Ref. |15] the SdH oscillations of a 2DEG confined in a gate-controlled InGaAs layer 
were observed under different gate voltages V g , at values -1V,-0.7V, -0.3V, 0V, 0.3V, 0.5V, 
and 1.5V. The electron effective mass ranges from m* = 0.049, at V^=-1V, to m* = 0.052, 
at V^=1.5V, and the corresponding sheet density changes from n — 1.6 to n — 2A1 x 10 12 
cm -2 . The first nodes B\ corresponding to the V g values given above are at fields Bi= 2.032T, 
2.025T, 2.011T, 1.98T, 1.923T, 1.894T, 1.87T, respectively, and the second nodes at B 2 = 
1.13T, 1.079T, 1.035T, 0.966T, 0.915T, 0.9T, 0.871T. Employing Eq. flUJ) we evaluate the 
Rashba parameter as a function of the gate voltage. We show the results in Fig. |U]as filled 
squares and fit them with the dashed curve, obtained with a = 7.04 — 2.26V^, + 0.87V^ 2 . Our 
results are consistent with those given in Ref. 3] obtained by fitting the first nodes (open 
circles) of the observed beating pattern with those obtained from an approximate evaluation 
of the resistivity; the fitting of the second nodes is shown by the filled circles. Our calculated 
magnetoresistivity, as a function of the magnetic field, for V g = 0.3V and T = 0.4K, is shown 
in the inset of Fig. |H1 Here we find the second node is well fitted and a smaller number of 
oscillations (n 2 — ni~ 26) than that observed (~ 28). It is worth noting that the value of a 
we obtained, after analyzing these two examples for InAs-based heterostructures, is of the 
same order of magnitude as that found by the microscopic model proposed in Ref. |l4| for 
comparable densities; from Fig. 3 of Ref. |l4j the extracted value of a is ~ 2 x 10~ n eVm 
at density n s = 10 12 cm -2 . 

Our way to extract the parameter a from the experimental SdH oscillations leads to 
theoretical results that are in rather good agreement with those obtained by fitting experi- 
mental curves. One advantage of Eq. (|2*4^) is that it is independent of the Zeeman splitting. 
Accordingly, the conclusion can be drawn that the Rashba effect plays the main role in the 
formation of beating patterns in the SdH oscillations in the measurements discussed above. 
However, as stated above some mismatches exist, e.g., the a value extracted from different 
node sets can vary and not all observed nodes can be fitted well with the same accuracy. 
Also, the measured dependence of the resistivity on the magnetic field is not well recovered 
by this simple model. This might be a result of the approximations introduced in it, e.g., 
the neglect of the bulk SOI in the model, the simplified impurity potential, the small u 
approximation or some unconsidered mechanism influencing the resistivity. 
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IV. CONCLUDING REMARKS 



We studied magnetotransport in a 2D electron system in the presence of the Rashba spin- 
orbit interaction term. When the subband broadening is much smaller than the Landau level 
separation, the effect of this term on the conductivity is manifested as a splitting of the SdH 
peaks. For week magnetic fields, with a level broadening comparable to the Landau level 
separation, a beating pattern appears in the conductivity plot as a function of the magnetic 
field. By measuring the position of two successive nodes, we can estimate the strength 
a of the Rashba term. The theory is in reasonably good agreement with the available 
experimental observations for p xx . In strong magnetic fields, where the integer quantum 
Hall effect is observed, a sufficiently strong a creates new plateaus between the integer 
plateaus in the Hall resistivity and splits the SdH peaks of p xx . 
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VI. APPENDIX 

The x and y components of the velocity operator read 



v x = 



dH 

dp x 



-ihV x /m* - u c y 



-ia/h 



ia/h 



-ihV x /m* - u c y 



(26) 



dH 

d P y 



—ihVy/m* a/h 



a/h -ihVy/m* 



(27) 



Setting £ s = uo c l c [s 1 ' 2 + V s V s+l {s + l) 1 / 2 ]/^, T s = u c l c [(s - l) 1 / 2 + V 8 V a . lS ^\/y/2, 

g s = lu c I c [Vs^s 1 ! 2 - V s (s + 1)V2 - V2(a/hio c l c )V s V s+1 ] /^/2A S A S+1 , and H s = 
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vJc [DgS 1 / 2 — V s _i(s — l) 1 / 2 — \j2al hu c l c ] /y/2A s A s -i, we can express the matrix elements 
of v x and v y in the Landau representation as follows: 



(K(k x )\v x \K,(k' x )) = - [£s ri ^ h] 6 s , s ,_ 1 5 kx , k , - l^^ W^ (28) 



s+1 



s-1 



, T /, m , T „,M « [£ a — aVg/h] r r i \F S — aV s _i /h] r 

(y;(k x )\v y \y;,(k' x )) = - [ a /nrnr ± S sS-iSk.M + 1 n^A^ % ,s>+i5k xK (29) 



s+l 



■\j J\l, S J\l. S — \ 



(*:(k x )\v x \*+(k' x )) = (K(k x )\v y \*:>(K))\^-c 



(30) 



(*t(k*)\vv\*$(K)) = (K(k x )\v y \*;,(k' x ))\ a ^ a 



(31) 



(32) 



The matrix elements of the position operator in the y direction are: 



(K(k x )\y\K'(K)) 



(33) 



t ,r j [* 1/2 + ( s + l) 1/2 )?W+i] V-i<Wi 



(*t(k x )\y\K'(K)) 



lr 



[V s V s+1 s 1/2 + (s + If' 2 ] 5 s , sl ^5 kxK 



\J2A s A s +i 



The form factors |F^/(-u)| 2 read 



\s'/s) l ' 2 LU(u) + P m P„,Lf- s ( M ) 



S! U 



e (36) 
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{s'/s) l l 2 V rn V rn ,LU{u) + L s s '- S (u) 



2 s! M s " 5 

5 ! +A.a*A.i 



(37) 
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